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[ Abstract. The Chandrasekhar limit for strange stars described by a linear equation of 
^ ■ state (describing quark matter with density-dependent quark masses) is evaluated. The 
■ maximum mass and radius of the star depend on the fundamental constants and on the 
O . energy density of the quark matter at zero pressure. By comparing the expression for 
the mass of the star with the limiting mass formula for a relativistic degenerate stellar 
. configuration one can obtain an estimate of the mass of the strange quark. 
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1. Introduction 

One of the most important characteristics of compact relativistic astrophysical objects 
is their maximum allowed mass. The maximum mass is crucial for distinguishing be- 
tween neutron stars and black holes in compact binaries and in determining the outcome 
of many astrophysical processes, including supernova collapse and the merger of binary 
neutron stars. The theoretical value of the maximum mass for white dwarfs and neutron 
stars was found by Chandrasekhar and Landau and is given by Afmax ^ (Jtj^b^^^^ 
( ^hapiro fc Teukolsky 1983 ), where ms is the mass of the baryons (in the case of white 



dwarfs, even pressure comes from electrons; most of the mass is in baryons). Thus, with 
the exception of composition-dependent numerical factors, the maximum mass of a de- 
generate star depends only on fundamental physical constants. The radius i?max of the 

f, / f, \i/2 

degenerate star obeys the condition i?niax < ( Gm^ ) ' '^ith m being the mass of 



either electron (white dwarfs) or neutron (neutron stars) ( ^hapiro fc Teukolsky 1985 ) 



White dwarfs are supported against of gravitational collapse by the degeneracy pressure 
of electrons whereas for neutron stars this pressure comes mainly from the nuclear force 



between nucleons (Shapiro & Teukolsky 1983). For non-rotating neutron stars with the 



2 Harko & Cheng: Maximum Mass and Radius for Linear EOS 

central pressure at their center tending to the hmiting value Pc(?i an upper bound of 



around ZMq has been found (Rhoades & Rufhni 1974). 

The quark structure of the nuclcons, suggested by quantum cromodynamics, opens 
the possibility of a hadron-quark phase transition at high densities and/or temper- 



atures, as suggested by Witten (Witten 1984). In the theories of strong interaction, 
quark bag models suppose that breaking of physical vacuum takes place inside hadrons. 
As a result, vacuum energy densities inside and outside a hadron become essen- 
tially different and the vacuum pressure on the bag wall equilibrates the pressure of 
quarks, thus stabilizing the system. If the hypothesis of the quark matter is true, 
then some neutron stars could actually be strange stars, built entirely of strange mat- 



ter ( Alcock, Farhi fc Olinto 1986 , Haenscl, Zdunik fc Schaeffcr 1986 ). For a review of 
strange star properties, see ( ICheng, Dai fc Lu 199^ ). 

Most of the investigations of quark star properties have been done within the frame- 
work of the so-called MIT bag model. Assuming that interactions of quarks and gluons 
are sufficiently small, neglecting quark masses and supposing that quarks are confined 
to the bag volume (in the case of a bare strange star, the boundary of the bag coincides 
with the stellar surface) , the energy density p(? and pressure p of a quark-giuon plasma 
at temperature T and chemical potential [if (the subscript / denotes the various quark 
flavors u, d, s etc.) are related, in the MIT bag model, by the equation of state (EOS) 
( pheng, Dai fc Lu 1998] ) 

(1) 

where B is the difference between the energy density of the perturbative and non- 
perturbative QCD vacuum (the bag constant). Equation (^ is essentially the equation 
of state of a gas of massless particles with corrections due to the QCD trace anomaly 
and perturbative interactions. These corrections are always negative, reducing the energy 



density at given temperature by about a factor of two (Farhi fc Jaffe 1984). For quark 
stars obeying the bag model equation of state (Q) the Chandrasekhar limit has been 
evaluated, from simple energy balance relations, in ( [Bannerjee, Ghosh fc Raha 200C ). In 
addition to fundamental constants, the maximum mass also depends on the bag constant. 

More sophisticated investigations of quark-gluon interactions have shown that Eq. 
(|l|) represents a limiting case of more general equations of state. For example, MIT bag 
models with massive strange quarks and lowest- order QCD interactions lead to some 
corrections terms in the equation of state of quark matter. Models incorporating restora- 
tion of chiral quark masses at high densities and giving absolutely stable strange matter 
can no longer be accurately described by using Eq. (|l|). On the other hand, in models in 
which quark interaction is described by an interquark potential originating from gluon 
exchange and by a density-dependent scalar potential which restores the chiral symmetry 



Harko & Cheng: Maximum Mass and Radius for Linear EOS 



at high densities (Day et aL 1998), the equation of state P — P [p) can be well approx- 



imated by a linear function in the energy density p ( Gondek-Rosinska et al. 200C). It is 



interesting to note that Frieman and Olinto (Frieman & Olinto 1989) and Haensel and 



Zdunik (Haensel & Zdunik 1989) have already mentioned the approximation of the EOS 



by a linear function (see also Prakash et al. (Prakash, Baron & Prakash 1990), Lattimer 



et al (Lattimer et al. 199C). Recently Zdunik (Zdunik 2000) has studied the linear ap- 
proximation of the equation of state, obtaining all parameters of the EOS as polynomial 
functions of strange quark mass, the QCD coupling constant and bag constant. The scal- 
ing relations have been applied to the determination of the maximum frequency of a 
particle in a stable circular orbit around strange stars. 

It is the purpose of this paper to obtain, by using a simple phenomenological ap- 
proach (which is thermodynamical in its essence), the maximum mass and radius (the 
Chandrasekhar limits) for strange stars obeying a linear equation of state. Of course the 
maximum mass of compact astrophysical objects is a consequence of General Relativity 
and not of the character of motion of matter constituents. However, the formulae for 
maximum mass and radius, due to their simple analytical form, give a better insight 
into the underlying physics of quark stars, also allowing us to obtain some results which 
cannot be obtained by numerical methods. For example, from the obtained relations one 
can find the scaling relations for the maximum mass and radius of strange stars in a 
natural way. 

The present paper is organized as follows. The maximum mass and radius of quark 
stars with a general linear equation of state is obtained in Section 2. In Section 3 we 
discuss our results conclude the paper. 



2. Maximum mass and radius for strange stars in the linear approximation 
of the EOS 



We assume that the strange star obeys an equation of state that can be obtained by 
interpolation with a linear function of density in the form: 



p = a(p-po)c^ 



(2) 



where a and are non-negative constants, po is the energy density at zero pressure. 
Such an equation of state has been proposed mainly to describe the strange matter 



built of u, d and s quarks ( Gondek-Rosinska et al. 200C, Zdunik 2000). The physical 
consistency of the model requires po > 0. 

The particle number density and the chemical potential corresponding to EOS (0) 



are given respectively by (Zdunik 2000) 
1 p \ 



n (p) = no 1 + 



a PqC^ 



(3) 
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m(p)=/.o 1 + ^^ , (4) 

where no is the particle number density at zero pressure and /io = p^(? jn^. 

The parameters a and po can be calculated, for realistic equations of state, by using a 



least squares fit method (Gondek-Rosinska et al. 200C, Zdunik 200C). For the equations 
of state incorporating restoration of chiral quark masses at high densities proposed in Dey 



et al. (Dey et al. 1998) one obtains the values a = 0.463, po = 1-15 x VS^'^gjcm?' and a = 
0.455, po — 1.33 X lO^^g/cm^, respectively ( pondek-Rosinska et al. 2000[ ). The standard 
bag model corresponds to a = 0.333 and po = 4 x X^^^gjcrv? ( pheng, Dai fc Lu 199S ). 

From Eqs. (^)-(|^) it follows that the particle number and chemical potential are 
related by the equation 

^ = • (5) 

ftp 

From the numerical studies of strange star models we know that the density profile 
of this type of astrophysical object is quite uniform ( plcndcnning 1996|) . Therefore we 
can approximate n k, N /V ^ which leads to 



/xo V 3 / V'^o 
where N is the total number of particles in a star of radius R and volume V . 

With the use of Eqs. (||)-(H) one obtains the energy density of the star in the form 

P - f^)""" + ^p„. (7) 

a+l\3/ \no J a + 1 

The total mass M of the star is defined according to M — An J^^ pr^dr w ^pB? and 
is given by 

po ^47^^-ViV^"+'._3a 47r a ^ . 



a + l\ 3/ \no J 3a+l 
where we assumed that the energy density is approximately constant inside the star. 

Extremizing the mass with respect to the radius R by means of dM/dR = gives 
the relation 



a+1 

, , o , , — ) R^""" = ^^—poR'- (9) 

a + 1 \ 6 J \no J 3a+l 

Substituting Eq. into Eq. (||) we obtain the maximum mass of the strange star in 
the linear approximation of the EOS: 
Air 

M=—poR'. (10) 

This expression is very similar to the expression for the maximum mass of the quark 
star obtained assuming that the star is composed of three-flavour masslesss quarks. 



confined in a large bag ( Bannerjee, Ghosh fc Raha 200C , Cheng fc Harko 200C ). From 
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a physical point of view, Eq. (|10|) describes a uniform density zero pressure stellar type 
configuration. 

The maximum equilibrium radius corresponds to a minimum total energy of the star 
(including the gravitational one), for any radius. For ordinary compact stars, the mass 
is entirely due to baryons, and the corresponding (Newtonian) gravitational potential 
energy is of the order Eg ^ —aGA'P/R {a = —3/5 for constant density Newtonian stars). 
For quark stars, assumed to be formed of massless quarks, the total mass can be calculated 
from the total (thermodynamic as well as confinmcnt) energy in the star. One possibility 
for the estimation of the gravitational energy per fermion is to define an effective quark 



mass incorporating all the energy contributions ( Bannerjee, Ghosh fc Raha 2000 ). 

The gravitational energy per particle (the strange star is assumed to be formed from 
fermions) is 

= (11) 

where m^ff is the effective mass of the particles inside the star, incorporating also 
effects such as quark confinment. For a star with particles one can write M = 
NrUeff = PaV 1 or ^e// = P^jn. On the other hand one can assume p — po/2n 



(Bannerjee, Ghosh & Raha 2000), leading to msff — 2p/c^ . Hence, with the use of Eqs. 
(|), @ and ^ we can express the gravitational energy per particle as 

Eg = -2(^X G^R^ (12) 



, 3 / N 

The energy density per particle of the fermions follows from Eq. (0) and is given by: 
_ 47r 1 po 3 

The total energy E per particle is 



Extremizing the total energy with respect to the radius (with the total particle number 
kept constant), (f§)^_co„gt — Oi it follows that the maximum radius of the equilibrium 
configuration in the linear approximation of the EOS is given by: 

Rmax = Ro . ■ (15) 

V^TT {a+l)Gpo 

The maximum mass of the star can be calculated from Eq. (|l^) and is: 

M -4 1 
*'"-^-3(„ + i)3/2Gv^- 

In Eqs. (|l|) and (|l|) i?o is a numerical factor of the order Rq w 0.474. 
The maximum radius of the quark star given by Eq. (^5|) is the radius corresponding 
to the maximum mass. On the other hand for the existing models of strange stars, 
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the configuration with maximum mass has a radius which is lower than the maximum 
radius. For example, for strange stars described by the bag model equation of state, the 
maximum radius is 11.40fcm, while the radius corresponding to the maximum mass is 
10.93/cm, which is 4% lower than the maximum radius. This difference is neglected in 
Eq. (|l|). 

The maximum mass and radius of the star are strongly dependent on the numerical 
value of the coefficient Rq and estimations based on other physical models could lead to 
different numerical estimates of the limiting values of the basic parameters of the static 
strange stars. 



3. Discussions and final remarks 

In the present paper we have shown that there is a maximum mass and radius (the 
Chandrasekhar limits) for quark stars whose equation of state can be approximated by 
a linear function of the density. We have also obtained the explicit expressions for Mmax 
and Rmax- 

With respect to the scaling of the parameter po of the form po — )■ kpo, the maximum 
mass and radius have the following scaling behaviors: 

For the maximum mass of the strange stars this scaling relation has also been found 
from the numerical study of the structure equations in the framework of the bag model 



( IWitten 1984 [Hacnscl, Zdunik fc Schaeffer 19861 ). 

A rescaling of the parameter a of the form a + 1 — s- ii" (a + 1), with po unsealed, leads 
to a transformation of the radius and mass of the form 

Rmax ^ ^ ^ " Rmax ^ max ^ ^ ^ ^max- (-^^) 

A simultaneous rescaling of both a and poi with a+l^if(a+l),po^ ^Po gives 

Rmax ^ ^ ^ ^ ^ Rmaxj -^^max ^ ^ ^ ^ ^ ^max- (-^^) 

The maximum mass and radius of strange stars with linear EOS is strongly dependent 
on the numerical value of po, the mass decreasing with increasing pQ. For po — AB, with 
the bag constant B = IQ^^g/cm^ {bGMeVfm^'^) we obtain Mmax — 1.83Af0, a value 
that must be compared to the value M^ax — SM© obtained by numerical methods 



(Witten 1984, Haensel, Zdunik & Schaeffer 1986). The difference between the numerical 
and theoretical predictions is around 10%. For po — 1.33 x lO^'^ gcm^^ the maximum 
mass of the star is about IMq. 



Generally our formulae (15) and ( |16D underestimate the maximum values of the mass 
and radius because we have assumed that the density inside the star is uniform. It 
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is obvious that near the surface the density is much lower than at the center of the 
compact object. Due to the approximations and simphfications used to derive the basic 
expressions, reflected mainly in the uncertainties in the exact value of the coefficient Rq, 
Eqs. ( p^ ) and ( |l6| ) cannot provide high precison numerical values of the maximum mass 
and radius for linear EOS stars, which must be obtained by numerically integrating the 
gravitational field equations. 

For the the linear EOS, Mmax and Rmax depend mainly on the fundamental constants 
c and G and on the zero pressure density po (the bag constant). The Chandrasekhar 
expressions for the same physical parameters involve two more fundamental constants, h 
and the mass of the electron or neutron. 

For quark stars, usually one assume they are composed of a three- flavour system of 
massless quarks, confined in a large bag. Hence the mass of the quark cannot play any 
role in the mass formula. But the linear EOS with arbitrary a can describe quark matter 
with non-zero quark masses (the mass of the strange quark rus ~ 200MeV), forming a 



degenerate Fermi gas (Gondek-Rosinska et al. 200C, Zdimik 2000). Therefore this system 
should also be described by the same formulae as white dwarfs or neutron stars, not only 
by Eqs. (|l|) - ([l|). Generally po is a function of the mass of the strange quark, so this 
mass implicitly appears in the expression of the maximum mass and radius. But on the 
other hand we can assume that the Chandrasekhar limit also applies to quark stars with 
the baryon mass substituted by an effective quark mass ruqef / , representing the minimum 
mass of the quark bubbles composing the star. Hence we must have 



Eq. (^0|) leads to the following expression of the effective mass of the "elementary" 
quark bubble: 

/fi 1/3 \'/^ 

m,eff - j . (21) 

The effective quark mass is determined only by elementary particle physics constants 
and is independent of G. From its construction mqef j should be relevant when the system 
is quantum mechanical and involves high velocities and energies. With respect to a scaling 
of the zero pressure density of the form po — > kpo, the effective quark mass has the 
scaling behavior mqejf niqeffk^^^, similar to the scaling of the strange quark mass 



( IZdunik 20001) . 

For Po = 4 X 10^'^ gcm,-^ we obtain m,e// ^ 3.63 x 10^^^ g w 20AMeV. For po 



1.33 X 10^^ gcm-^ Eq. (21) gives m,e// ^ 4.9 x lO^^s^ 27bMeV. The mass given 



by Eq. (21) can be considered as the minimum mass of the stable quark bubble. It is 
of the same order of magnitude as the mass of the strange quark. Therefore the 
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Chandrasekhar limit applies also for quark stars if we take Wqe// for the mass of the 
elementary constituent of the star . 

In the present paper we have considered the maximum mass and radius of strange 
stars in the linear approximation of the equation of state and the dependence of these 
quantities on the parameter a has been found. We have also pointed out the existence of 
scaling relations for the maximum radius of strange stars, an aspect that has not been 



mentioned in previous investigations (Witten 1984, Haensel, Zdunik fc Schaeffer 198(:, 



Bannerjee, Ghosh fc Raha 2000| , [Zdunik 20001) . Our formulae also lead to the transfor- 
mation relations for the maximum mass and radius of strange stars with respect to 
separate and simultaneous scaling of the parameters a and pq. On the other hand the 
possibility of estimation of the mass of the strange quark from general astrophysical 
considerations can perhaps give a better understanding of the deep connection between 
micro- and macro-physics. 
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